ACCEPTED BY ApJ 

Preprint typeset using IAT^X style emulateapj v. 20/04/00 



O 

(N 

< 



O 

u 

43 

Oh- 
i 

o 

CO 



> 

(N 

cn 
O 

o 



13 



A DOUBLED DOUBLE HOTSPOT IN J 0816+5003 AND THE LOGARITHMIC SLOPE OF THE 

LENSING POTENTIAL 

Katherine M. Blundell 1 , Paul L. Schechter 2 , N. D. Morgan 2 , Matt J. Jarvis 3 , Steve 

Rawlings 1 , & John L. Tonry 4 

accepted by ApJ 

ABSTRACT 

We present an analysis of observations of the doubly-lensed double hotspot in the giant radio galaxy 
J 0816+5003 from MERLIN, MDM, WIYN, WHT, UKIRT and the VLA. The images of the two hotspot 
components span a factor of two in radius on one side of the lensing galaxy at impact parameters of 
less than 500 pc. Hence we measure the slope of the lensing potential over a large range in radius, made 
possible by significant improvement in the accuracy of registration of the radio and optical frame and 
higher resolution imaging data than previously available. We also infer the lens and source redshifts to 
be 0.332 and > 1 respectively. Purely on the basis of lens modelling, and independently of stellar velocity 
dispersion measurements, we find the potential to be very close to isothermal. 

Subject headings: gravitational lensing: strong radio continuum: galaxies galaxies: quasars: individual 
(J0816+5003) 



1. INTRODUCTION 

In contrast to multiply-imaged, compact flat-spectrum 
sources, strongly-lensed extended, steep spectrum sources 
occupy a distinctly inferior position in the gravitational 
lens pantheon. Yet, in a paper that marks the beginning 
of the modern era in gravitational lensing, I Press &: Gunnl 
(|l973h argued that radio lobes were the most promising 
targets in searching for lensing. Indeed, examples in their 
figure 6 resemble the lens which is the study of this paper, 
but since that paper was written, lensed lobes have proven 
more difficult to interpret. 

Chief among the difficulties is the continuing poor reg- 
istration between the radio and optical reference frames 
(jDeutschl Il999() . Radio lobes tend not to have optical 
emission that allows their position to be determined 
with respect to any lensing galaxy that might be iden- 
tified. As lobes are extended and polarized, they offer 
more potential constraints than multiply imaged point 
sources, but modelling such emission requires substantial 
mach in ery dKochanek. Blandford. Lawrence fc Naravanl 
19891: lEllithorne et all |199fit iDve fc Warrenl 120051: 
Vegetti fc Koopmansl I2009T T For the most part, lensed 



lobes have been passed over in favor of their compact 
flat-spectrum cousins. 

We describe in this paper Merlin and VLA observations 
of the doubly-lensed double radio hotspot J 0816+5003 
(a.k.a . 6C 0812+501) discovered both by iLehar et al.l 
(|2001h . using a clever automated strategy for finding ex- 
tended radio emission which is gravitationally lensed, and 
independently by the Oxford group searching for high red- 
shift radio galaxies. The techniques for high redshift radio 
galaxy searches include filtering on t he basis of steep radio 
spectra and small angular size (e.g . iBlundell et allll99Rt 
iJarvis et al.ll200lt iJarvis et"lill200ll ) in the course filtering 
radio source surveys to optimize seaches for high redshift 



radio galaxies. J 0816+5003 is an example of a candidate 
included in the initial search but subsequently excluded by 
the filtering criteria. 

The lens we study in this paper is the hotspot com- 
plex within the northern lobe of the giant radio galaxy 
J0816+5003, which has an angular e xtent of over an ar - 
cminute and is pictured in figure 1 of ILehar et ail (l200l . 
Both the associated radio core, a half arcminute away to 
the south west, and its optical counterpart are identified, 
allowing accurate registration of the lensed hotspot and 
the lensing galaxy. In addition we find, wit h higher res- 
olutio n radio imaging than was available to ILehar et al.l 
(2001J), that the hotspot has two components both of which 
are doubly imaged. Finally, the lensed images are asym- 
metrically located with respect to the lensing galaxy, per- 
mitting determination of the radial profile of the lensing 
potential. In particular, the images of the two sources 
span a factor of two in radius on one side of the galaxy at 
impact parameters of less than 500 pc. 

This, together with accurate radio and optical astro- 
metric registration as well as higher reso l ution radio 
images than were available to ILehar et al.l ( 20011 ). en- 
ables us to make a pure-lensing measurement of the 
slope of the potential of the total mass of the lensing 
galaxy on the scales probed, including the dark mat- 
ter. By combining lensing and dynamical analyses, the 
SLA CS team have investigated this and other properties 
(e.g. IBolton et al.ll20q7l: iGavazzi et alJl2007t iBolton et al.l 
120081: iBarnabe et al.l |2009T >. The slope of the poten- 
tial figures prominently both in measurements of the 
Hubble constant and i n discussions of the mass profiles 
of elliptical galaxie s (IChen. Kochanek. fc Hewitt! 1995: 



Romanowskv & Kochanek 



Bernstein fc Fischerl 119991: 

1999 t IWilliams fc Sahal 1 2 00ft IWitt et all 120001 : ICohn et a] 
20011 ). There are very few examples where one can 
measure the slopes in this way - - some exceptions 
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include MG 1654+ 13 4 (iKochaneld I1995D. 095 7+56 1 
(|Grogin fc Narayanl Il996allbl: iFischer et all Il997ft . 
MG 2016+112 (iTreu & Koopmansl I2002D. B 1152+199 
(|Rusin et al ll2002fl. but cf. B 1 933+503 (ICohn et alll200lD 
and especially PG 1115+080 (IChen. Kochanek. fc Hewittl 
ll995tlRomanowskv fc Kochanekl 19991) . Some newly found 
lenses by the SLACS team have also yielded constraints 
on the slope, albeit in most cases by a comparision 
of stellar velocity dis persions and Einstein ring radii 
(jKoopmans et al.l(2006l) . 

With regard to measurement of Hq, a very small frac- 
tion of systems are suitable for measurement of time delays 
but as mentioned above, only a small fraction of systems 
are suitable for measurement of the slope of the poten- 
tial. Since the overlap of these two sets is slim, it remains 
important to glean as much as possible about the slopes 
of galaxies potentials and extrapolate fr om the former to 
the latter ([Kochanek fc Schechterl I2004D . The currently 
accepted value of H is reconcilable with inferred values 
from the 10 firm time-delay measurements in gravitational 
lenses only if dark matter halos are eliminated and a con- 
stant mass-to-light ratio used. 
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Fig. 1. — VLA contours of the north hotspot structure of 
J 0816+5003 at 5 GHz. It has a beam of 0.55" X 0.36" at a p.a. 
of 82.39 deg. The lowest contour is at 0.4mjy/beam and adjacent 
contours differ by a factor of y/2. The greyscale is an i?,-band image 
fro m the WHT . The convention for naming follows that instigated 
bv lLehar et~aTI J2001T) . 



2. OBSERVATIONS 

2.1. VLA imaging 

J 0816+5003 was observed for 2 minutes on 1996 Dec 06 
in the A-configuration of the Very Large Array (hereafter 
VLA) at 5 GHz under program code AR365 originally as 
part of an Oxford survey for high redshift radio galaxies. 
The data were reduced using standard procedures in the 
AIPS software package. The primary flux calibrator was 
3C 286 and the phase calibrator was 0804+499. The image 
showing its lensed northern hotspot structure is shown in 
Figured] 



2.2. Merlin imaging 

To obtain higher resolution radio imaging than in Fig- 
ure!]] J 0816+5003 was observed with six antennas of 
the Multi-Element Linked Radio INterferometer (here- 
after MERLIN) at 1.7 GHz on 1997 Jun 27 for 10.58 hrs, 
during which observations of a nearby phase calibrator, 
B 0803+514, were interleaved. Again, the data were re- 
duced using standard procedures within AIPS including 
standard polarimetric calibration made possible by the ex- 
tensive parallactic angle coverage of the phase calibrator 
throughout the 10+ hour observation. 

Figures^] and |3] show contours of the total intensity 
MERLIN data in different ways as described in the fig- 
ure captions, with the highest resolution image having a 
point spread function with FWHM of 0.08" x 0.08". 



2.3. Optical and near-IR imaging 

In order to establish the nature of any galaxy associated 
with 6C 0816+5003, this source was observed for 5 minutes 
in i?-band on 1997 Jan 9 on the William Herschel Tele- 
scope (hereafter WHT) using the Aux Port Camera, and 
was flat-fielded and reduced using standard procedures. A 
greyscale of this image is presented in FigureQ] 

Near-infra red K-band observations of J 0816+5003 were 
made using the UFTI instrument for 27 minutes on the 
United Kingdom Infra-Rcd Telescope (hereafter UKIRT) 
on 2002 Jan 13. Unfortunately conditions were not photo- 
metric, but approximate flux density boot strapping was 
made via a short (9-minute) service observation on 2002 
Mar 11. For each of these observations, following stan- 
dard practice, each successive 60-sec exposure was jittered 
from the previous position and a flat-field was created by 
combining all frames from a particular observation using a 
median filter . All data reduction was performed in IRAF 
and followed IWillott et ail (<2003D . 
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Fig. 2. — MERLIN image at 1.7 GHz of part of the lensed north 
hotspot in J 0816+5003 showing the southern two components of the 
three depicted in Fig[T] This image was made using the IMAGR rou- 
tine of the AIPS software package with a Briggs' robustness parame- 
ter of —3 and super-resolved with a restoring beam of 0.08" X 0.08". 
The lowest contour is 0.3mjy/beam and adjacent contours differ by 
a factor of \/2. The beam is 0.08" X 0.08". 
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Fig. 3. — MERLIN image at 1.7 GHz of the same part of the 
north hotspot shown in Fig[2] but with the image compressed in 
the East-West direction by a factor of four. The lowest contour is 
0.3mjy/beam and adjacent contours differ by a factor of y/2. The 
beam is 0.08" x 0.08". 

2.4. Redshifts of lens and lensed galaxy 

The ISIS spectrograph on the WHT was used on 2002 
Feb 5 using the red and blue arms of the ISIS spectrograph 
(with R158R and B300B gratings respectively), in order 
to find the redshifts of the background source, namely the 
radio galaxy, and the lensing galaxy. The spectra were 
reduced using standard IRAF procedures. The redshift of 
the lensing galaxy is inferred to be 0.332, on the basis of 
an observed 4000A break and the Ca H and K absorption 
lines (at 3968 and 3933 A respectively) shortward of this. 
At this redshift, 1" is 4.738 kpc assuming Hq = 71km s~ , 
r2o = 1 and A = 0.73. The background source, being 
a radio galaxy, belongs to a relatively homogeneous cat- 
egory of sources with respect to its host galaxy. This is 
displayed most app arently by the tightne ss of the radio 
galaxy K-z relation ()Lillv fc Longair|[l984lL I ndeed, with 
work on the radio galaxy K-z Hubble relation (jJarvis et all 
[200lL IWiliott et all[200l we are able to make a crude es- 
timate of the redshift of the radio source. A magnitude of 
K ~ 17.8 from the UKIRT observations implies a redshift 
of z>1.0, which is also consistent with its broad-band op- 
tical colours. We can also gain more information from the 
optical WHT spectrum, where the presence of a relatively 
bright continuum but the absence of any of the bright for- 
bidden narrow-emission lines, implies that the radio galaxy 
probably lies within the redshift desert at 1.2 < z < 1.7 
(i.e. [Oil] 3727 has been redshifted out of, and Ly— a has 
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yet to enter, the optical passband). While the if-band 
magnitude indicates that the radio galaxy is probably at 
z>l, the WHT spectrum tentatively indicates that z ~ 1. 
Thus, we can be certain that the radio galaxy is at a much 
greater distance than the lens. 
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Fig. 4. — Spectrum from the ISIS spectrograph on the WHT, 
observed on 2002 Feb 05. The redshift of the lens galaxy is inferred 
to be 0.332, on the basis of the H and K absorption lines shortward 
of the 4000 A break. 

3. OPTICAL ASTROMETRY 

To deduce accurately the slope of the potential of the 
lensing galaxy requires secure relative astrometry between 
the images and the lenses. In the case of PG 1115+080, 
a difference in the image positions of several parts in a 
thousand makes a difference between a best-fit power-law 
slope of a = — 1 rather than a = (jSchechteij [2000) , 
which relates to the surface mass density profile with ra- 
dius r as E(r) oc r a ~ 2 . In the case of J 0816+5003, the 
lensed images are only detectable in the radio waveband 
and the lensing galaxies is not visible in the radio. Thus, 
making deductions about the slope of the power-law of the 
galaxy lensing its hotspot depends on how well the radio 
and the optical/near- IR frames can be tied together. The 
derivation of astrometric solutions using an image from a 
telescope with reimaging optics or with a prime focus cor- 
rector (such as the WHT, UKIRT or Keck) and with a 
reference from the USNO catalog would not be sufficiently 
successful, due to off-axis distortions. 

A 17' x 17' i?-band image of this field was observed with 
the Michigan-Dartmouth-MIT 1.3 m telescope (hereafter 
MDM-1.3m) by Jules Halpern on 2002 Feb 22. The as- 
trometry of this image was fitted with John Thorstensen's 
automated astrometric program^ which matched 156 
bright stars over this field from the USNO catalog and 
put the astrometric solution into the fits headers. This 
astrometric solution was transferred to a 4096 x 2048 pixel 
i?-band image obtained from the WIYN with 0.5" seeing 
on 2002 Mar 14. 



The WIYN and MDM scales were both determined from 
fits which assumed that the scale was the same in both di- 
rections. The scale factor in both cases was determined to 
better than 1 part in 10 4 . The postional rms between the 
WIYN and MDM objects was 0.34 pixels which at a scale 
of 0.1406189"/pix is roughly 0.05". Thirty six stars were 
used in the WIYN fit. The error in the galaxy position 
was taken to be 0.05 pixel (as judged from the scatter in 
the different position measurements) that translates into 
a 0.007" error in the lens position (in each direction). The 
Right Ascensions and Declinations of the lensing galaxy 
and of the radio galaxy core and their relative separation 
were measured using an average of three differences ob- 
tained from DAOFIND, IMEXAM and VISTA and are 
reproduced in TableQ] 

The positions of the pair (E and W) of doubled (Nl 
and N2) hotspot components are listed in the first four 
rows of Table[T] and were determined using JMFIT within 
AIPS on the image shown in Fig[3] The positional uncer- 
tainty along each axis is 4 milli-arcsec for the unresolved 
N2-E, N2-W and Nl-W components, and because of be- 
ing slightly resolved that for the Nl-W component is 10 
milli-arcsec. The east and north offsets are measured with 
respect to the position of the radio nucleus of the radio 
galaxy (which is at 08 16 37.40774 +50 03 38.9782) which 
locks the radio and optical frames together because of the 
coincidence of the radio nucleus with the centre of the op- 
tical host galaxy. The positional uncertainty of the lensing 
galaxy is good to 0.007" and dominated by the 0.5'/ seeing 
of the WIYN image. The flux densities at 1.7 GHz of the 
compact lensed components are also measured using JM- 
FIT and are for the brightest compact features, excluding 
surrounding smooth extended emission. 

By applying a correction to the Right Ascension and 
Declination of the radio galaxy to align it with the radio 
core, and applying this exact same correction to the for- 
mer RA and Dec exact relative offsets between the lensed 
images and lens were then input to the model. 



4. DIMENSIONLESS MODEL FITTING 

The most northerly feature of the hotspot complex (N3) 
shown in Figure [1] is likely to be distorted (tangentially 
sheared) by the presence of the lensing galaxy; it exhibits 
no multiple compact components in the MERLIN data. 
Thus, as is often the case, the hotspot in this giant ra- 
dio galaxy h as multiple intrinsic (unlensed) components 
(|Laing| [l989). However, its impact parameter from the 



galaxy means that it is not multiply imaged (see Fig[S]). 

5 http: / /www. astro. lsa.umich.edu/obs/mdm/manuals/Manuals/automatch. html 
6 http: / /ftp. nofs.navy.mil/projects / pmm/ 



Doubly lensed hotspots and the logarithmic slope of the potential 



5 



component 



N2-E 
N2-W 
Nl-E 
Nl-W 
lens 



RA 
J2000 

08 16 38.82360 
08 16 38.81384 
08 16 38.86700 
08 16 38.82578 
08 16 38.82580 



Dec 
J2000 

+50 04 08.1275 
+50 04 08.2098 
+50 04 06.4285 
+50 04 06.2481 
+50 04 07.7950 



east offset north offset flux density 
arcsec arcsec mJy 



13.5881 
13.4944 
14.0046 
13.6090 
13.6412 



29.1493 
29.2316 
27.4503 
27.2699 
28.8530 



Table 1 

Fitted positions and flux densities of radio components. 



1.12 ± 0.11 
0.93 ± 0.11 
4.01 ± 0.12 
5.74 ± 0.12 
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Fig. 5. — Upper panel: the source plane showing as blue crosses 
the positions of the intrinsic double hotspot feature. Lower panel: 
image plane. 

4.1. General Approach 

Gravitational lenses can be modelled in terms of a di- 
mensionless two-dimensional potential <f> that is a function 
only of position on the sky. Dimensionful quantities - dis- 

7 http: / /redfive. rutgers.edu/~keeton/gravlens/index. html 



tances and masses - can then be derived from redshift 
measurements and specification of the Hubble constant. 

While there is in principle a great deal of information to 
be extracted from multiple images of extended sources, 
this quickly leads into dividing the sou rce into pixels, 
each o f which is taken as a f r ee pa rameter (|Dve fc Warrerj 
120051: IVeeetti fc Koopmansl I2009D . Radio hotspots are 
particular challenging bec ause there i s relatively little 
symmetry to be exploi t ed (lLainsJ [l989t iBlack et all 119921 : 
iHardcastle et al. 1 11997k iLeahv et allll997l ). Moreover im- 
age infidelities may be present in the details of the bright- 
ness structure. For example, while the resolution of the 
MERLIN image is approaching that ideally needed for this 
system, the six antennas of the MERLIN array undersam- 
ple much of the extended structure seen in the VLA im- 
age, whose resolution is insufficient to reveal the under- 
lying structure. Resolution of this awaits the emerging 
generation of radio telescopes especially the e-MERLIN 
and EVLA. For the present we limit ourselves to mod- 
elling the positions and fluxes of the compact components 
of the hotspots. 

Reasonable and robust first-order results can be ob- 
tained using only point sources. We use the excellent 
public domain program lensmodel written and gr aciously 
made available by Chuck Keeton (Kccton 2001)0 Our 
sources are clearly extended, but we try to take out the 
tangential stretching, seen clearly in the MERLIN image 
shown in Figure^ by squashing by a factor of four in the 
East- West direction. The resulting image can be seen in 
Figured Obtaining accurate radio positions of the double 
images of each pair were made easier by the fact that one 
of the knots shows no signal in Stokes Q at the resolution 
of the MERLIN data, so the position for the other pair was 
measured from a Stokes Q image. To measure the other 
pair, the Stokes Q image was subtracted from a polarised 
intensity map (y Q 2 + U 2 ). This approach overcame the 
problem that the image shown in Figure[5] does not com- 
pletely resolve the two blobs in N2. [N.B., the rotation 
measure in the lensing galaxy implied by the polarisation 
angle variation with frequenc y is only 20 rad/m 2 , modest 
by normal galaxy standards ( Leahy et allll986| ).] 



4.2. Details of the model 

We chose to model the two dimensional projected grav- 
itational potential with the function 

a/2 



Q/2 

y 



(1-ef 



(1) 
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where 9' x and 9' are angular position on the sky and e is 
the ellipticity of the equipotentials. In the limit of zero 
cllipticity, the radius of the Einstein ring, te can be ob- 
tained from the parameter b with te = 6 1 / 2 ~ Q . Potentials 
due to circularly symmetric isothermal galaxies (i.e. with 
a = 1) with external tides yield similar results but with 
slightly inferior fits to the data. For the small ellipticities 
considered here there is no advantage to using elliptical 
isodensity contours, which in the limit of small e are a 
factor of 3 flatter than the equipotentials. 



350 



300 



150 



100 



. Positional data \ 

2 

. Flux density \ 
.total x 



1.0 

a 



Fig. 6. — Plot of x 2 as a function of a, the power-law slope of the 
potential. 

The results of stepping through a and allowing all the 
other 11 parameters vary are shown in Figure[6] Our best 
fit is for a = 1.031, with an uncertainty, as judged from 
the change in x 2 , of 0.002. The best fitting model param- 
eters are given in TableEJ Note, in particular, that the 
best fitting position for the lens differs from our measured 
position by less than 0.003". 

5. SLOPE OF LENSING POTENTIAL 

The use of a particular gravitational lens system in de- 
ducing Hq is only as good as our knowledge of the lensing 
potential. The galaxy responsible for lensing the hotspot 
in J 0816+5003 has a power-law which is remarkably close 
to isothermal. 

In addition, it is unlikely that the lensing potential is 
much shallower than isothermal because ther e is no hint 
of ev i dence for any odd (e.g. 'third') image (W inn et al.l 
120031: iRusin k Mall2001ft . The flux density ratio of the 
peak of the brighter northern (N2-E) blob on the MER- 
L IN image to 3 x t he rms noise on this image exceeds 10. 

iKochanekl (| 19951 ). similarly, finds that in MG 1654+134 



that 0.9 < a < 1.1. IRusin et all (12002ft find in B 1152+199 
that 0.95 < a < 1.21 and iTreu fc Koopmani (|2002t ) find 
in MG 2016+112 that a = 1.0 + 0.1 + 0.1. However, seem- 
ingly inconsistently with t h is em erging general picture, 
IChen. Kochanek. k Hewittl (|1995f) claim that the mass 
profile of the lensing galaxy in MG 1131+0456 to be in- 
compatible with a = 1. I n add ition, for PG 1115+080 
Ro manowskv fc Kochanekl (|l999t ) report that no viable 
lens model can yet be ruled out. 

6. FUTURE WORK AND CONCLUSIONS 

Sensitive radio images at resolution better than a few 
10s of milli-arcsec would considerably broaden the scope 
of a lensing model to be fitted, and this is within the capa- 
bility of current technology. The key observational limita- 
tion however, remains the registation of the optical frame 
with respect to the radio frame and widespread efforts to 
improve this situation generally should remain a priority. 

We conclude by remarking that our pure-lensing mea- 
surement of slope of the potential of the lensing galaxy 
in this system as being very close to isothermal is a mea- 
surement of the total mass on the scales probed, includ- 
ing the dark matter. However, it sheds no light on the 
"bulge-halo conspiracy" that separately the dark and light 
matter do not appear to be isothermal — the light matter 
being steeper than isothermal, the da r k matter being shal- 
lower dTreul 120071; iBolton et all [20071 120081: iGavazzi et ail 
120071: iBarnabe et al.ll2009T) — but collectively they appear 
isothermal over an increasing range in redshift as the sys- 
tem we have studied in this paper illustrates. 
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measured e from n. in the limit of zero ellipticity the radius of the elnstein ring, in arcseconds, is given by 

r=6 l/(2-a) < 
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